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Planetary systems based on a quantum-like model 
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Planetary systems have their origin in the gravitational collapse of a cloud of gas and dust. 
Through a process of accretion, is formed a massive star and a disk of planetesimals orbiting the star. 
Using a formalism analogous to quantum mechanics (quantum-like model), the star-planetesimal 
system is described and the flow quantizing the gravitational field theoretical model parameters 
are obtained. Goodness of fit (chi-square) of the observed data with model quantum-like, to the 
solar system, satellites, exoplanets and protoplanetary disk around HL Tauri is determined. Shows 
that the radius, eccentricity, energy, angular momentum and orbital inclination of planetary objects 
formed take discrete values depending only on the mass star. 


I. INTRODUCTION 


In the solar system the distance of the planets from the Sun follow a simpl e ge ometric pro gressio n, known as 
the Titius-Bode law: r = 0.4 + 0.3 x 2 n AU with n = —oo, 0,1,..., 7 Lynch 20031 iNeslusan 2004| . However, 


better agreement using the Bohr equation is obtained: r = 0 .0425 x n 2 AU, with n = 3,4,5,6,8,11,15,21,27 
[Caswell 19291 iPenniston 1930L IBarnothv 1946al IBarnothv 1946bj . Comparisons of the Titius-Bode law with the 
equation of Bohr raised the issue of the applicability of certain principles of quantum mechanical in systems or¬ 
bital of the astronomy [Corliss 1986l |. To interpret Boh r’s equatio n within the framework of classical mechan- 


to the Jovian satellites 

Rubcic & Rubcic 1995, Rubcic & Rubcic 19981 

Hermann et al 19981 and extrasolar planets 

Rubcic & Rubcic 1996, 

Rubcic & Rubcic 19991. Finallv. it has attempted to describe planetary systems with the 

Schrodineer equation Reinisch 1998, de Oliveira Neto et al 20041 Smarandache & Christianto 2006, Nie 2011|. 


In quantum mechanics, Planck’s constant is extremely small, h ~ 10~ 34 J.s and the ratio de Broglie wavelength is 
inversely proportional to momentum. For a massive object length de Broglie wave is very small, and consequently the 
wave behavior of a macroscopic object is undetectable, making quantum mechanics reduces to classical mechanics. 
Thus quantum mechanics is a theory developed to explain phenomena only a microscopic scale. 

Planetary systems originate from a nebula which collapses to form a star and a disk of planetesimals orbiting the 
star. The trajectory of the planetesimals is determined by the principle of stationary action, SS = 0. However, 
planetesimal is disturbed by the presence of the other, oscillating around the classical trayectory; this causes for a 
closed and stable trayectory, the action is quantized; that is, the action is an integer multiple of a quantum of action, 
S = nh s (see Appendix A). The quantum of action, h s (which plays a role equivalent to Planck’s constant h) is a free 
parameter to be determined and depends on the physical system in question. Considering the flux quantization of the 
gravitational field, the value of h s is determined. For a massive object: h s h, this allows to describe a macroscopic 
system with the formalism of quantum mechanics (quantum-like model). This means that the formalism of quantum 
mechanics is applicable to any system that quantize the action. One could use the Bohmian interpretation to give 
meaning to the formalism. 

This implies that planetary systems (planets-Sun, satellites-planet, exoplanets-star, exoplanets-pulsar) must quan¬ 
tize not only the orbital radius as demonstrated, but also eccentricity, orbital inclination, angular momentum and 
energy. Through the goodness of fit (chi-square) shows that there is a very good agreement between the observed 
data and the quantum-like model, although the formulation is based only on the interaction between two particles. 


II. FLUX QUANTIZATION OF THE GRAVITATIONAL FIELD 

In a microscopic system, the electric potential energy between a proton and an electron, is given by, U e = — e 2 /47re 0 , 
and an equivalent macroscopic system: the gravitational potential energy between a star and a planetesimal, is given 
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by, U g = —GM s m/r. A relationship of proportionality between the electric potential energy and the quantum of 
action h (Planck’s constant), is assumed: U g ~ h, also between the gravitational potential energy and h s parameter: 
U g ~ h s , which plays the same role of Planck’s constant, in a macroscopic system. We can establish the relationship 
between these two energies, to obtain a dimensionless independent amount of r, and normalize the macroscopic 
parameter on a unit, the Planck constant: U g /U e ~ h s /K = H. 

The flow of the gravitational field generated by the mass of a star, through a spherical gaussian surface, is given 
by: <hg = 47 tGM s and the electric field flux generated by a proton is: 4> e = e/e 0 . If we take the relation between the 
intensity of the gravitational and electromagnetic force in terms of their flows, $ = mQ g /e$> e , we have, $ ~ H. It is 
considered <f>, as flow of the quantum of action per unit solid angle: $ = H/Q (II acts as the proportionality constant), 
we obtain 


-lnGM s fi h s 
e 2 /4ireo m 

the Bohr radius is given by: a s = (h s /m) 2 /GM S , and the energy of the ground state: E s = GM s /2a s . We can express 
these quantities in terms of the solar mass: H s = kh.Q , a s = ka@, and E s = kEQ, where k = M s /Mq called the 
scale factor. The theoretical parameters are given by: ft 0 /m = 7.62311 x 10 14 Js/kg, u 0 = 2.92705 x 10~ 2 AU and 
Eq/tti = —15.1538 GJ/kg. 


III. ORBITAL DYNAMICS MODEL 


In classical mechanics, the trajectory followed by a planetesimal, orbiting a star, is an ellipse with semi-major axis a 
and eccentricity e. The average value of the distance planetesimal-star during a complete orbit, turns out to be (r) = 
a(l + e 2 /2). The total mechanical energy is E/m = —GM s /2a and the magnitude of the orbital angular momentum 
is L/m = y 7 GM s a( 1 — e 2 ). By conservation of angular momentum, in the accretion process of planetesimals, planets 
must be in the same plane 0 = 0. The equation of motion for the ith planet orbiting a star is: f\ = — G(M s +nii)/r 2 +Fi , 
where Fj are small magnitudes, which contain the perturbative effects of all other objects on the ith planet, which 
is not a point particle; this causes the parameters which determine the orbit, vary periodically. The variation of the 
orbital kinetic energy, due to the variation of the eccentricity is: A E^i f'j/m = GM s /2 (1 — € 2 )/rf — (1 — e 2 )/?y 
and, due to the variation of the orbital inclination: A Eg^jGm = (2GM s /ri) sin 2 [(0/ — 0j)/2]. 

In the quantum-like model, by substituting the Hamiltonian, H = —h 2 X7 2 /2m — GM s m/r (V 2 is the Laplacian 
in spherical coordinates) in the Schrodinger-like equation, we obtain: E n /m = —E s /n 2 , where E s = GM s /2a s 
and a s = (H s /m) 2 /GM S . Also, the average distance, (r n ,t) = a s /2 [3n 2 — £ (1 + 1)]; the magnitude of the angular 

momentum, Li/m = sjt (£ + l)h s , and its orientation, 0£ im( , = arccos (ni // \Jl (( + 1) j. The quantum numbers 
leading to orbit and where the probability is maximum, are n = 1,2,..., £ = n — 1 and me = ±£. The sign of me 
corresponds to normal and retrograde orbits, respectively. With the wave function: 


V'n = 


n 2 y/{2 n — 1)! 


fir 1 (M) 


(where p = r/a s ), the probability of finding a planetesimal in space | pi\) n (p, 9, <j >)| 2 is obtained (see Figured]). 

If the planetesimal obeys the quantum-like model, depending on the distribution of matter and the process of accretion 
of planetesimals, a large object (planet, dwarf planet or asteroid) is formed on the mean value of the distance 
(rt) = a s n(n + 1/2), ie, the object must orbit the star formed, following an elliptical trayectory with semi-major axis 

at = a s n 2 , eccentricity e t = 1 / y/n and inclination 0 t = arccos \J\ — 1/n (because this is the most probable, and stable 


closed orbit); with a total energy, E t /m = — E s /n 2 and magnitude of the angular momentum, Lt/m = \/n(n — 1)A<,- 
The orbital inclination and eccentricity of the formed objects changes cyclically over time Lisiecki 20101 ]. so the 
observed total energy E 0 /m = —GM s /2a 0 should be corrected: 


E c — E 0 + A E e ( t ^ — A + Ae. 


(1) 


in this equation Ae corresponds to other effects, which generally lead to a loss of energy, such as planet-asteroids 
collisions and others. 
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Figure 1: Probability of finding a planetesimal at a point in space, \ptpn (p, 9, </>)| 2 , for (a) n = 1 and (b) n > 1. 


IV. THE SOLAR SYSTEM 

On the origins of the solar system, the formation of objects by the process of accretion of planetesimals can occur 
anywhere in the protoplanetary disk, but the stable closed orbits give rise to larger objects, according to the quantum¬ 
like model, this occurs in at = kaQn 2 . Based on the uncertainty principle, ApAat — h s , we have that A at = kaQ. For 
an average radius of ~ 300 km, the icy moons and rocky asteroids in our solar system passing from a rounded potato to 
a sphere [Lineweaver 20ld |. for this reason, we selected objects with a diameter greater than Vesta > 525.4 km (NASA 
JPL Small-Body Database), smaller objects are discarded because they can easily change orbit, due to collisions or 
internal dynamic processes (eg, Yarkovsky effect). Assuming that the scaling factor is k = 1, assigning a quantum 
number to each object and selecting the of greater diameter: (4) Mercury, (5) Venus, (6) Earth, (7) Mars, (9) 
Vesta, (10) Ceres, (13) Jupiter, (18) Saturn, (26) Uranus, (32) Neptune, (36) 78799 (2002 XW93), (37) Pluto, (38) 
Haumea, (39) (2010 KZ39), (40) Makemake, (41) (2013 FZ27), (42) 42301 (2001 UR163), (44) 84522 (2002 TC302), 
(45) (2013 FY27), (46) (2010 RE64), (48) Eris, (50) 229762 (2007 UK126), (56) 145451 (2005 RM43), (58) (2008 
ST291), (95) (2012 VP113), (134) Sedna, (Figure O- The filling of the orbits depends mainly on the mass of the 
star and the amount of material available in the protoplanetary disk. The Bohr radius value observed is obtained: 
a.Q = (2.92007 ± 0.00256) x 10 -2 AU, Xd 0 f = 0.30574, R 2 = 0.99997, and percentage error e% = 0.24 with respect 
to the theoretical value. 

Figure [3] shows the relationship between the corrected energy E c , equation ([TJ) , and the quantum numbers. Once 
have formed interior objects orbital radius is maintained, while objects in the energy continuum (> 30 AU) can be 
easily moved (Kuiper belt). Jupiter is at the inflection point. Is obtained E^/m = (14.65188 ± 0.00160) GJ/kg, 
XdoF = 0.00002, R 2 = 0.99956, and e% = 3.31. It is assumed that the percentage error is due mainly to the effect 
they have had the impact of asteroids on the inner solar system objects. In turn, the degree of fit, is evidence of 
the quantization of angular momentum, ie, the quantization of the orbital inclination and magnitude of the angular 
momentum in the period of formation of the solar system. 

In Figure U shows the observed angular momentum, which depends on the eccentricity, L 0 /m = yjGM s a 0 (1 — e„); 
the objects far away from the Sun have very high orbital inclinations, this generates a variation in the orbital 
eccentricity m, which in turn produces an anomalous behavior in the angular momentum. Suppressing effect of 
orbital inclination on the eccentricity, we can obtain a corrected angular momentum: L c /m = yjGM s r 0 (1 - e 2 ), is 
obtained: h° Q /m = (0.76144 ± 0.00081) PJs/kg, Xn„F = 0.03682, R 2 = 0.99992, and e% = 0.11. 

The Sun corresponds the quantum number n = 1, consequently, its angular momentum is zero and the probability 
distribution of matter is a sphere, see Figure[IJa), this could explain why the angular momentum of the Sun corresponds 
only ~ 2% of the entire solar system. The quantum numbers n = 2,3,..., have a probability distribution of matter 
in the form of concentric toroids, see Figure [1(b), which have different degrees of inclination and orbital velocity, 
the absence of planets in numbers quantum n = 2, 3, indicates that objects migrated outwards or planetesimals were 
captured by the Sun. 
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Figure 2: Orbital radius (planets): theoretical, at (solid line) and observed, a 0 (circles). 



Quantum number, n 


Figure 3: Energy (planets): theoretical Et (solid line), observed E 0 (gray circles), and corrected E c (black circles). 


A. Orbital migration 

The parameters obtained correspond to the current solar system, but the planets had a different distribution when 
they formed: it is assumed that the sequence of the inner planets was, (4) Mer cur y, (5) Venu s, (6) Earth, (7) Mars, 
(8) Vesta and (9) Ceres. The Jovian planets has the sequence given by (Gomes et al 2005j . In this solar system 
objects representing at least 26 is considered that at the time of its formation, the same objects were present. Taking 
the maximum goodness of fit (maximum likelihood) for these conditions we obtain the sequence: (4) Mercury, (5) 
Venus, (6) Earth, (7) Mars, (8) Vesta, (9) Ceres, (13) Jupiter, (16) Saturn, (19) Neptune (21) Uranus, (34) 78799 
(2002 XW93), (35) Pluto, (36) Haumea, (37) Makemake, (38) 55565 (2002 AW197), (39) (2010 RF43), (40) 42301 
(2001 UR163), (41) 84522 (2002 TC302), (42) (2004 XR190), (43) (2013 FY27), (44) (2010 RE64), (45) 225088 (2007 
OR10), (46) Eris, (48) 229762 (2007 UK126), (53) 145451 (2005 RM43), (55) (2008 ST291), (90) (2012 VP113), (127) 
Sedna and a° Q = (3.25085 ± 0.00279) x 10~ 2 AU, x 2 Do f = 0-29805 R 2 = 0.99997, and a scal e factor k = 1.1106 2, 
therefore the mass of the sun was M s ~ 1.1 M 0 , which agrees with [Boothrovd et al 1991 Ij . and [Guzik fe Cox 1 995], 











5 



Figure 4: Orbital angular momentum (planets): theoretical Lt (solid line), observed L 0 (gray circles), and corrected L c (black 
circles). 



Figure 5: Orbital radius (satellites): theoretical at (solid line) and observed a 0 (black circles). 


B. Satellites 

Is have selected the satellites of the planets of the solar system (Planetary Satellite Physical Parameters, JPL, 
NASA) with diameters > 350 km. In order to compare with the solar system, has taken the ratio of the orbital radius 
and the scaling factor ( a 0 /k ). By the same procedure, used for the solar system, is obtained: (10:J) Io, (12:S) Mimas, 
(13:J) Europe, (14:S) Enceladus, (15:S) Tethys, (16:J) Ganymede, (17:S) Dione, (21:S) Rhea, (21:J) Callisto, (23:N) 
Proteus, (31:S) Titan, (32U) Ariel (37:U) Umbriel, (48:U) Titania, (53:S) Iapetus, (55:11) Oberon, (171:E) Moon. Was 
excluded (40:N) Triton because anomalous behavior with respect to the model, shows that is a captured by Neptune 
object. The values ??obtained for the parameters are observed: a q = (2.92645±0.00173) x 10~~ 2 AU ( XdoF = 0-26477, 
R 2 = 0.99999, and e% = 0.02 (see Figure [5]). Is obtained: Tig/m = (0.76581 ± 0.00185) PJs/kg, Xn oF = 0.14702, 
R 2 = 0.99982, and e% = 0.46; E° /m = (15.11796 ± 0.02155) GJ/kg, xjjoF = 1 - 2Q8 x 10 " 7 : r2 = 0.99993, and 
e% = 0.24. 
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Quantum number, n 


Figure 6: Orbital radius (exoplanets): theoretical at (solid line) and observed a 0 (black circles). 



Figure 7: Number of exoplanets vs. difference between the orbital radius theoretical at, and observed a 0 . 


V. EXOPLANETS 

T his researc h has made using the Exoplanet Orbit Database and the Exoplanet Data Explorer at exoplanets.org 
[Wright et al 201 1| . Have been selected systems with a single star with known mass and exoplanets, radio and orbital 
eccentricity known; assigning a quantum number and selecting the most massive, there are 219 exoplanets. In order 
to compare with the solar system, has taken the ratio of the orbital radius and the scaling factor (a 0 /fc). The value 
obtained for the Bohr radius is: Sg = (3.03542 ± 0.00944) x 10 -2 AU, Xd 0 f = 0.00108, R 2 = 0.99777, and e% = 3.70 
(see Figure 0). 

In the Figure [T] is shown the number of exoplanets vs. difference between the ob serv ed and theoretical orbital 
radius (A a/k), has made an adjustment to a Gaussian distribution y = y 0 + (. A/w\J it/ 2) exp [—2(A a/kw) 2 ], here 

y 0 = 2.66240 ± 0.88034, A = 10.01864 ± 0.20650, w = 0.04750 ± 0.00094, a/k = 0.00517 ± 0.00072, X d 0 f = 2.32095, 
and R 2 = 0.99978. This means that exoplanets are not formed anywhere, only preferentially in the orbits given by 
at/k = aQn 2 . 

In classical mechanics, by conservation of angular momentum, the objects must be in the same plane, 0—^0° (in the 
accretion disk), in the quantum-like model there is a quantization of angular momentum and its orientation, objects 
have higher orbital inclination between the closer you are to the star, this explains the elevated orbital inclinations 
found for exoplanets (0 —> 90°). It have /Iq/to = (0.76619 ± 0.002910 PJs/kg, Xd oF = 0.00681, R 2 = 0.99638, and 
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Figure 8: Orbital angular momentum (exoplanets): theoretical Lt (solid line), observed L 0 (gray circles), and corrected L c 
(black circles). 


e% = 0.51. 

The total energy shown large variations when the exoplanet is very close to the star (before continuous), this 
is due to variations in the eccentricity (AE e ( t ,o)) and the inclination (A Eg( t ,o)) of the orbit, is obtained Eq/vi = 
(17,57860 ± 0.36720) GJ/kg, X doF = 0.24013, R 2 = 0.88301, and e% = 16.00. 


A. Protoplanetary disk around HL Tauri 


Considering a centra l protostellar mass for HL Tau given by k — 0.55 Beckwith et al 1 990. 


ISargent, fc Beckwith 199il |. has to be formed planets in at = fca©n 2 . Based on uncalibrated image of the proto¬ 
planetary disk HL Tauri taken by ALMA (ESO / NAOJ / NRAO) where planets may be forming, several gaps are 
distinguished a t = (10.06, 24.49, 46.94, 68.02) AU which correspond to n = (25, 39, 54, 65), respectively. 


VI. CONCLUSIONS 

The quantum-like model is applicable to any system that quantize the action, therefore, the formalism of quantum 
mechanics can be used in macroscopic systems where, the quantum of action, h s h. In this particular case, this 
occurs because the flow of gravitational field is quantized. The dynamics is determined only by the mass of the central 
object (as the case, star or planet) independently of the mass of the orbiting object. The probability distribution 
of planetesimals in space and the physical quantities (radius, eccentricity, energy, orbital angular momentum and its 
inclination) of the objects in the planetary system, take a certain value (quantized). Therefore, planets and exoplanets 
tend to be formed at a predetermined distance which depends on the mass of the star and the quantity of matter 
on the disc, obeying at = kaQti , the orbit may be normal or retrograde. Similarly, the angular momentum and its 
orientation takes certain value, however, because the planets and exoplanets are not punctual objects, occur periodic 
variations, in the eccentricity and inclination orbital, causing a change in the value of the observed energy, equation 
©. The obtained theoretical model can be used as a tool to parameterize and study the formation of planetary 
systems. 
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Appendix A: QUANTUM-LIKE MODEL 

In the protoplanetary disk particles condense into small planetesimals of mass m, which orbit the star. Since the 
mass of the star, M s , is extremely large compared to the mass of planetesimals, M s to, action star prevails and 
the interaction between them, only results in weak perturbations to their respective orbits. The Hamiltonian (in 
polar coordinates) for a planetesimal orbiting a star is: H = p 2 /2m + U(r), where p r = mr and U(r) is an effective 

potential energy U(r ) = L 2 /2mr 2 — GM s m/r, in turn L 0 = pg = mr 2 9 is the angular momentum. As the force is 

conservative, F = —S7U', planetesimal experience a attractive force, for r > r a , and repulsive force, for r < r a , while 
r = ro forces cancel each other out. It can easily be shown that U(r ) has a minimum at r = r 0 : U 0 = U(r) r — ro 
necessarily (VU) r=ro = 0 and k 0 = (V 2 t/) r=r > 0. 

The lower energy state corresponds to E = U 0 , planetesimal path is a circumference of radius r 0 centered on 
the star. The planetesimal on an orbital cycle runs its perimeter ( x 0 = 2m' 0 ) at a time ( t 0 ). As, Vt/ 0 = 0 the 
planetesimal in this kind of orbit equivalently behaves a free particle, in this case the equation of Hamilton-Jacobi 
is given by: (VS') 2 /2m + dS/dt = 0. Is denoted by x the perimeter distance traveled by the planetesimal. As in 
the Hamiltonian coordinate x (cyclic), or the time t does not appear explicitly, there are two conserved quantities: 
the linear momentum (VS = p 0 ), and the total energy of the system (— dS/dt = E 0 ). Consequently the action is 
separable: S(x, t) = S x (x) + St(t) = p 0 x — E 0 t, and is periodic: S(x,t) = S(x + x 0 , t + t 0 )■ 

Bertrand’s theorem shows that under any initial condition, the only potential producing stable orbits have the 
functional form ~ 1/r (Kepler) or ~ r 2 (harmonic oscillator); additionally closes the orbit if the relationship between 
the radial and orbital frequency is w/w 0 = n/m, where n and to relatively prime. By disturbing the orbit of a 
planetesimal it can be shown that the planetesimal oscillates around its original path. A perturbed orbit if after n 
periods of variation of r (t a = nr), the complete revolutions planetesimal m, nA8 = to27t is repeated and a revolution 
in the planetesimal travels a distance A = A 9r 0 , where Bohr’s rule is obtained for the orbits, nX = 27r r 0 . 

Hamilton’s principal function vanishes for: 5(0,0) = S{x 0 ,t 0 ) = n(p 0 A) — n(E 0 r) = 0. The terms inside the 
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parentheses correspond to a constant which is denoted by h s . Note that the action is quantized, that is, S x (x 0 ) = St(t 0 ) 
is given by n times a quantum of action h s , we will call macroscopic parameter, which must be determined for this 
physical system. The quantum of action h s when the planetesimal travels a x a = X is h s = S x ( A), or when time 
elapses t a = t is, h s = 5*(r); these expressions are analogous to the relations of de Broglie and Planck, respectively: 


p 0 = h s k, E a = h s uj, (Al) 

where, h s = h s /2n called, reduced macroscopic parameter, k = 2 tt/X and to = 2 it/t. 

The perturbation of the trajectory of a planetesimal is caused by the presence of other planetesimals, causing it to 
move from its original path r D , a distance £ = r — r 0 (which is very small, r — r 0 ). Making a Taylor series expansion of 
the effective potential U ( r) around of r Q , its radial motion is given by: p^/2m + fc c £ 2 /2 = E — U 0 , where p^ = ml ; is a 
transverse momentum with respect to the unperturbed motion. Consequently, the disturbance makes the planetesimal 
harmonically oscillate around the non-perturbed trajectory. 

Based in Hamilton canonical equation, p £ = —dH/d£, we may obtain the wave equation, d 2 ip(x,t)/dt 2 — 
v 2 X7 2 ip(x,t) = 0; where ip(x,t) = £,{x)t;(t) is defined and v 0 = ix/k. The corresponding solution has the form: 
ip = C exp [i(k x — ui t)\. As the path is periodic, ^>(0,0) = ip(x 0 ,t 0 ) = 0 then: k —> k n = rnr/X, ui —> ui n = titt/t 
being (n = 1,2,...). As we can see, the boundary conditions have the effect of discretizing the wave number and 
frequency. The solution of the wave equation takes the form: ip n {x,t) = C n exp [i{k n x — oj n t)\. 

However, substituting relations m in ip n (x,t) there is no correspondence between the wave speed and the speed 
of planetesimals. Planetesimals to represent using the wave equation, it is necessary to modulate the amplitude of 
the wave with a positive real function, A(x, t) = A(x — v 0 t), to limit the extent of the wave and cause the envelope to 
move with a group velocity, v Q . Obtaining, in general: 


’k (r, t) = C A(r, t) exp 



(A2) 


the exponential allows taking into account the phenomena related to the superposition of waves, while the coefficient 
corresponds to a distribution function (or probability density) to find the planetesimal in space: g = |fk(r, t)\ . By 
integrating over all space the complex constant C can normalize to unity. The probability current is defined as 
J = gv 0 . Conservation of probability is given by the continuity equation: 


+ v J = o, ( A3 ) 

the probability density moves in space with the same velocity X/S/m = v 0 and follow the path set of planetesimals. 
The Hamilton-Jacobi equation we add the term Uq = \k 0 £, 2 which corresponds to the effect of the disturbance, 

H + UQ + ~dt =0 ' (A4) 

this implies that the classical force is being affected by another force, which generates the quantization of the action 
of the physical system: F = —X7U — X7Uq. The eikonal equation: V 2 A + ( 2mUQ/h 2 )A = 0 allows us to define a 
relationship between Uq and A. Substituting this expression in equation (IA4I) and with (IA3I) we can construct an 
analogous equation to Schrodinger, which call Schrodinger-like: 


H^ (r, t) = ih a (r, t) . (A5) 

ot 

From the Hamilton-Jacobi equation (IA4I) we obtain the trajectory followed by the particles and from Schrodinger 
equation (IA5I) we obtain the distribution of particles in the space or the probability density, p(r, t). 

In the star-planetesimal system planetesimals that quantized action giving rise to stationary waves (or stable 
trajectories) which store the energy of the disturbance. These waves interfere constructively giving rise to phenomena 
of resonance or stationary states (when the system with the resonance frequency is disturbed), are described by 
the eigenvalue equation Hip(r) = E n p(r). Planetesimals not quantized action are represented by traveling waves 
transporting the disturbance energy, these waves end transferring its energy, so that the planetisimals assume the 
classical trajectory, which is described by (1A4I) where Uq = 0, ie, the formalism of classical mechanics. Therefore, we 
have a privileged system where states appear when Uq ^ 0, which can be described with the formalism of quantum 
mechanics. 




